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1. INTRODUCTION 

The final report on our LANDS AT-2 studies in Switzerland has to “be reviewed 
in connection with our ERTS-1 and S3CYIAB-EREP final reports (lO.l).All these . 
studies were carried out continuously by more or less the same group on in- 
vestigators and using the same methodological approach* 

Various topics were studied from LANDS AT- 2 (and additional) data. Main emphasis 
were concentrated on snow surveys and on the development of a semi- operational 
system to automatically map and measure the areal extent of the snow cover and 
its changes based on digital MSS-data (7*l). 

In addition land use studies were undertaken in the Po Valley around Milano 
and in the Swiss Plateau, here with special reference to applications for 
planning purposes* Of particular interest are the studies on the correlation 
between areal extent of settlements and population (7*2)* 

The investigation will not be terminated with this final report. Various sub- 
stantial projects (including four dissertations and two H.A. thesis) are under 
progress. The National Science Foundation of Switzerland just agreed to support 
the studies again for another three years. Therefore the activities will con- 
tinue on 'an even broader scale dealing with quite different topics (8). 

The main research, primarely on digital snow mapping and on data processing 
problems in general, is undertaken in a joint effort by the Department of 
Photography, Swiss Federal Institute of Technology (Prof .Dr. W.F, Berg) and the 
Department of Geography, University of Zurich (Prof .Dr.H.Haefner) , funded ’oy 
the Swiss National Science Foundation. 

In particular the following persons were or are engaged in this research: 

- Seidel, Klaus, Ph.B., project leader 

- Itten, Klaus, Ph.D., snow mapping, land use planning, system analysis 

- Binzegger, Robert, Ph.D., land use mapping 

- Lichtenegger, Jiirg, M.A., land use mapping, mirror experiment 

- Stanz, Karl , M.A., atmospheric effects, ground measurements 

- Muri, Ren 6, M.A., snow mapping, evaluation of various classification algorithms 

- Stirnemann, Hans-Peter, snow mapping 

- Fasler, Fritz, Ph.B., system specialist (part time) 
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- Funkhouser, Arthur, M.A., control functions of photomat ion system 

- Matt, Otto, system specialist (part time) 

- Paschke, Anton, system specialist (part time) 

In addition various other institutions were actively engaged in the LANDSAT-2 

project, namely: 

- Swiss Institute for Snow and Avalanche Research, Weissfluhjoch-Davos 
(j. Martinec, Fn.D.) 

- Department of Geography, Swiss Federal Institute of Technology, Zurich 
(Prof .D.Steiner) 

- Department of Geography, University of Berne (prof «B •Messer li/K.Winiger, Ph*D.) 

On the other hand we have to point out some problems, which hampered the progress 

of our research vrork considerably and which should be borne in mind when reading 

this final report: 

- Since it vras never known to us, if our test area was covered during a specific 
orbit, it was very difficult to arrange for connecting underflights and field 
observations. Nevertheless during part of the melting season in 1975 and most 
of 1976 ground observations were carried out at each overflight and aerial 
coverage was taken when the sky was clear. 

- Data arrived always with a lengthly timel&ps, which unabled a direct compa- 
rison of the image information with the ground conditions (especially regar- 
ding position of snow line). 

- Additional difficulties arose from the switch in data delivery from GSFC to . 

/ 

TELESPAZIO, Rome. During spring, summer and fall of 1976 we did not receive 
any data from the Pacino station. The first tapes arrived around November 8th, 
1976. This was a handicap to our research work, in particular to the mirror 
experiment (7.2.2), since we could not check whether the mirror orientations 
were correctly executed. 

- In total, there are only a few cloudfree scenes of our snow mapping test 
site in Grisons (Pig. 2) during, the melting periods, which do not provide 
a sufficient sequence to systematically map the changes of the snow cover. 

Only the latest receipt of data from TELESPAZIO as well as GSFC provided 
us a more or less continuous successtion of interpretable scenes for the 
melting period 1976 (Fig. l). They are being processed right now. 



- 5 - 


The situation for land use classification in the Swiss Plateau (test site 
"Grosses Moos", Kt. of Berne, Fig. 2) is somewhat better covered for the 
vegetation period of 1976, but again we got these data only very recently. 

The author expresses his appreciation to Dr.K.Itten, Dr.K. Seidel, J.Lichten- 
egger, M.A. , K.Stanz, K.A., R.Muri, M.A. , and Miss E. Leuzinger for their 
helpful assistance in preparing the report. 



2* OBJECTIVES 


The longterm objectives are to develop an integrated classification system to 
operationally man and measure dynamic surface features . in particular snow 
and land use, by means of multitemporal digital satellite data. Applications 
of snow mapping should lead to better surface water runoff predictions and 
provide a better data base for hydrologic balance. Satellite land use surveys 
will make available the most recent situation for planning purposes. 

To monitor change detection in an automated procedure it is inevitable to 
include geometric corrections for a direct combination of the data from 
different orbits. Radiometric corrections are needed to compare data, taken 
under various atmospheric conditions , being from the same orbit when classi- 
fying large areas (e.g. on the North or South side of the Alps) or on multi- 
temporal scenes. Since no other studies deal with this problem under European 
conditions, which differ substantially from the ones in North America* research 
on these problems was introduced. 

Studies on geometrical problems deal also with the merging of LANDS AT data with 
data from other origin such as from maps (e«g. with a digital terrain model or 
a geographical information system) or ground measurement. As a first step into 
this direction we set up an experiment to exactly locate ground points on the 
LANDSAT data matrix by marking specific pixels with mirrors. 

Methodological research concentrates on data processing , e.g. the modification 
of, existing software for our purposes and the development of nev/ programs. In 
the last months it was possible to install an interactive image analysis system 
PDF ll/40 with RAKTBK-display-imit at the Swiss Federal Institute of Technology 
in Zurich* which will be an integral part of our data processing system, it is 
hoped that the interactive approach will especially improve, facilitate and 
accelerate the selection and testing of the training samples for a supervised 
classification procedure. 

Besides it was our constant concern to not only develop an adequate classifi- 
cation system but also to dispose of a high quality cutout device . This could 
be achieved by procuring a photomation system P 1700 by OPTRONICS Inc., e.g. 
an electro-optical drum scanner and writer. With this system it is not only 
possible to present the results of the classification on b&w film in a pixel 
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by pixel representation but also to create color composites* A special study 
is under way to investigate all electronic and optical transfer functions of 
the photomation system, to set up extensive control programs, and to enlarge 
the software programs* 

A secondary field of application for LANDS AT-image s are studies in developing , 
countries , which are carried out at the Department of Geography University 
of Zurich* Objectives as well as the methodological approach are quite diffe- 
rent from the digital snow mapping project. They shall only be summarized 
here briefly* 

The objectives are to interpret natural geographic features, in particular 
geological and geomorphological surface elements, the hydrological network 
and the vegetation cover* The interpretation is done on the imagery by con- 
ventional methods using b&w as well as especially produced color composites 
to improve the interpretability. 

The first study on a “qualitative geological-geomorphological interpretation 

p 

of Sfil’S-l multispectral images" covered. 110' 000 km in the Hoggar-Hassif , 
Sahara, and evaluated the possibilities of LAKDSAT-data under arid conditions, 

A second study deals with the regionalization of the Yemen Arab Republic, based 
on maps of the landform features, the surface stream network (delineation of 
watersheds etc.), and the vegetation pattern, which will improve substantially 
the present knowledge on this developing country. She research was undertaken 
in form of two K.A. thesis. 
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3* SUMMARY OP DATA OBTAINED 

Since the start of LANDS AT— 2 7 data of Switzerland was obtained as summarized 
in Pig. I. To achieve a complete coverage of Switzerland, two frames each 
on three following days are necessary. Up to now no total coverage of Switzer- 
land on three following orbits was covered in one sequence. 

Since March 1976 our test site — the Swiss Alps - could he covered again by 
LANDSAT-l, received by the Fucino station. At the same time the activities 
of G-SPC continued with LANDSAT-2 over Europe to a certain extent. Therefore 
it was possible to get some time sequences of nine days. But again, due to 
technical reasons as well as weather conditions, only a very limited number 
of such data exists (Pig. l). 

The quality of the data is sufficient in general; regarding the digital tapes, 
which we just received from TELESPAZIO, it is very satisfactory. 
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4, GROUND TROTH AlfiD AERIAL UNDEJRPLIGHTS 


4 * i Snow mapping 

To support the digital snow mapping project, extensive measurements and 
field observations were carried out during each LANDSAT-2 overflight at 
various points within our tost site in Central Grisons (Fig. 2). In addition 
high-altitude aerial underf lights could be arranged on cloudfree days with 
the Swiss Air Reconnaissance* An example of such a high-altitude b&w aerial 
photography is given in Fig, 5, 

On the ground the following parameters were noticed or measured: 

- position of temporary snow line (visible observations, photography, mapping) 

- type of snow 

- characterization of snow surface (forms, cleanness, wetness, age etc.) 

- characterization of snow pack (granularity, wetness, density etc.) 

- snow temperature (at surface) 

- snow depth 

- snow mass 

- air temperature 

- air moisture 

- slope angle and exposure 

- cloudiness 

The purpose of these ground investigations is a twofolded one* First the 

• X 

information is used for a better delineation of the test samples for the 
supervised classification. In this connection we studied the problem of the 
most important snow parameters, ¥hich ones of the above mentioned parameters 
have a specific influence on the spectral signature, and consequently should 
be measured furthermore as a minimum requirement ? As a result of this field 
campagne we learned that for a description - of the physical properties of 
the snow or the snow -surface all above-mentioned parameters are needed. 

Only air temperature and density are highly correlated. 

Second the informations, in particular on the exact location and course of 
the transidient snow line, are used to verify the accuracy of our snow classi- 
fications. 

Measurements of the spectral properties of various snow types were also 
carried out. The results of this stud;/ are summarized in 7.1.1. 
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4.2 Land use change detection (mirror experiment) 

Between May 15th and October 30th 1976, 15 field surveys were undertaken 
in the test site "Grosses Moos 1 ’ (Big. 2) during LAEfDSAT— 1 and —2 overpasses. 
The following duties were conducted: 

- Setting and adjustment of four mirrors with theodolite and alidade (7.2.2) 

- Sampling of surface conditions, soils, land use, phenology etc.; 

- EXOTECH-measurements of the most important land use categories (corn, 
summer— wheat , winter— wheat , rye, meadows). The instrument was mounted 
on a long ladder about 3? 5 meters above the vegetation canopy. 

In addition a detailed land use map was compiled of the entire test area. 
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5. HARDWARE 

Over the last four years the already existing hardware at the Swiss Federal 
Institute of Technology, Zurich, and at the University of Zurich, could be 
substantially enlarged and improved for image analysis and data processing 

purposes . 

♦ 

Today the following hardware is available: 

- IBM 570/155 

- CDC 6400/6500 

- Interactive image analysis system PDF II /40 with RAMT3K display-unit (not 
fully operational yet) 

- Mini-Computer HP 9850 

- Photomation system OPTRONICS P-1700 

- Kicrodensiiometer system with automated data registering 

- Quantimet QTK 720 

- Digitizer D-HAC 

- Radiometer EXOTECH-lOO 

- Stereoplotter KERN PG-2 

- Zoom Transfer Scope ZT4-A 



6 . METHODOLOGICAL APPROACH 


There are three aspects to be considered in digital classification of multi- 
. yactral scanner data: 

- preprocessing 

- feature extraction 

- presentation of the classification results 

Fig. 4 gives a theoretical approach to decision making techniques . It is based 
on the fact that the spectral signature of each ground element represents a 
four-dimensional "features space", where each element is characterized by a 
single point. In order to classify pixels the feature space must be subdevided 
into different ,r sub-vo lumes" or categories, either by a supervised technique 
based on training samples or by a non-supervise d one, applying clustering. ' 
The different decision making techniques are summarized in Fig. 4 * 

6.1 Preprocessing 

Two different tasks have to be considered here. The first one deals, with the 
reformating and organization of the digital MSS data for our own hardware com- 
ponents and our specific purposes. The second one contains all corrections of 
the data in accordance to certain criteria to improve the classification. 

These programs include geometric as vzell as radioinetric corrections and re- 
arrange the data for an output in a maplike form. 

Since Dec. 1975 all parties' interested' in digital image processing around 
Zurich decided on a new data format (Zurcher Datenformat, ZDF) to coordinate 
and facilitate the communication and the interchangeability of the data bet- 
ween the various hardware systems and user groups. 
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6*2 Feature extraction 

Up to now we have concentrated on a supervised approach (Fig 1 . 4), based on 
training samples. Since the accuracy of the classification depends heavily 
on the correct choice of these training samples; they are delineated and 
statistically tested with specific care. The training samples are actually 
subject to ©tcvfcistioal fluctuations. Thb l&rger the number of pixel# bad 
the closer together they are, the “better their definition. Using the 
HP-9830 with a plotting device it i 3 possible to receive a two-dimensional 
graphical presentation of the distribution of the samples in the feature 
space. Combining always two different variables it is possible to get a good 
idea of the position and the separability of the training groups within the 
feature space (Fig. 5 )* 

To improve the discrimination, new synthetic variables such as the r&tioing, 
were introduced, in particular for land use interpretation ( 7 . 2 . 1 ). To avoid 
to include pixels which are too far away from any one training sample to make 
a sensible assignment, thresholding is used and the corresponding pixels 
classified as unassignable* Future steps for a simpler and more automated 
procedure of sample selection are directed toward the installation of a 
11 spectral signature bank 11 and the use of the interactive system (2 ). 

Various algorithms were tested such as 

- D-Class (distance in feature space) 

A-Class (angular position in feature space) 

- Maximum likelihood 

- Stepwise linear discriminant analysis 
^ P PD (parallel epiped method) etc. 

The usefulness of some of these algorithms (accuracy versus' expenditure) was 
tested with FREP-S-192 data for snovr mapping purposes (MURX, 1976) and will 
be investigated in more detail in the near future for various geographical 
applications. 

For the IBM an own classification, system based on a PPD-algorithxn was developed 
which allows a fast and easy interpretation, even by untrained people. Both, 
an interactive (IBIS) system and one for batch mode (BIS), are existing 
(Pig. 6). 

An additional system evaluation of several well-known interactive image pro- 
cessing systems was undertaken for snow mapping purposes (7.1.3) . 
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6.3 Presentation of results 

* 

One of the great problems in digital processing is an adequate presentation of 
the classification results as images or even in maplike form, Present-day 
standard computers are illequipped for this type of output. Photos taken from 
TV monitors often lack the necessary geometrical accuracy. Special equipment 
is required -for high-quality presentations. 

The available instrument at Zurich is an electro-optical drum scanner and 
writer (Photomation System P 1700 by OPTRONICS Inc.), which records pixel 
by pixel on a sheet of film (Fig. 15). Prom these b&w films it is easy to 
make color composites. . 

During the reporting period a software system has been created (which is sub- 
ject to further improvement) which is directed 

- to a more simple and user-oriented manipulation of the system and 

- toward a greater flexibility in the data-output products. 

Today it is possible with the photomation system to 

- select any part out of a frame 

- reproduce this part in a linear or non-linear fashion 

- influence the grayscale reproduction (MTF) 

- create a grid system on the image 

- analyze the image data 

- produce histograms of a magnetic -tape file (e,g. for an optimal tone 
reproduction) 

- add synthetic information on the output film for controlling and identifi- 
cation purposes, such as step wedge (sw), registration and location marks 

« i 

(KR), text in alphanumeric form (LT). 

Furthermore feature extraction is possible with a PPD-algorithm, whereby the 
class boundaries can be introduced directly for each channel. This procedure 
has great economical advantages for ^classifying large areas, because it 
allows direct feature extraStioa and output in pictorial form without having 
to use a big computer, : 

Future activities will include a careful evaluation of the optical transfer 

i ft 1 ; 1 * . 

functions of the system and the integration of more sophisticated classifi- 
cation algorithms based on the t ab le- 1 o ok-uo-aopr o a ch . 
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7. ACCOMPLISHMENTS AND PROBLEMS 

Two broad areas of applications are studied: snow mapping for runoff predictions 
and land use for planning purposes. 


7.1 Snow mapping 

By far the most important factors for snowmelt runoff prediction by a day to 
day calculation or using mathematical models are always the extent and the 
changes of the snowcovered areas. It is the aim to not only monitor the 
changes of the snow cover hut also to gain evidences on the snowpack and the 
coniditions of the snow surface. For these purposes ground measurements on the - 
spectral behavious of snow, were conducted. 


7.1.1 Spectral reflectance of snow 

During winter and melting period 1975/76 measurements of spectral properties 
of snow were conducted on the ground and from a helicopter with an EX0TECH-100. 
It was found to be more appropriate to measure under various meteorological 

conidtions in the field than under controlled lab situation. Measurements of 

r \ ) 1 

the reflected radiances were made of snow surfaces at different slope angles 
in sun and in shadow of six different snow types/ Measurements from a heli- 
copter served as means for investigating atmospheric radiation attenuation 
of the lower atmosphere. Both types of measurements wiil he used for a corre- 
lation with IANDSAT brightness values. N % 

Up to now the results show that two main groups of snow are presumably discer- 
nible using the four MSS bhnds (Fi’g* 7). She most important snow parameters 
(as observed in the field and described in STAENZ, 1976), modifying the re-, 
fleeted radiances^ are: grain form;- snow depth (up to 50 cm); snow density, 
wetness, pollution of surface, roughness of surface, age of: snowpack. 

h * ' • . r i. 

Pig 8 and 9 give examples -.of the target reflectance in. function of snow depth 

\ 1 ' 

for channel 4- and channel ,7 respectively. There is a clear distinction bet- 
•' I r . r ; 

ween fresh, moist, snow; old, wet but clean snow; and old, wet and polluted 

M '* f \ 

snow. .The influenpe of the snow wetness on the brightness values was measured 

.. • >• 'I •? u * , 

during about two ^ours at ;the same place. The snow surface changed from frozen 

U i \ 

to wet during the measuring period, From the target reflectance values measured 
in all four channels, ratios were compiled.^ They are presented in Fig, 10 in 
function of time. 
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Ratio 4/7 and to a lower amount 5/7 show a distinct modification from the frozen 
and dry to the wet snow conditions. The combination of these two channels con- 
sequently will provide good informations on the moisture conditions of the snow 
surface. A comparison of the reflectance measurements in sun and shadow are 
given in Fig. 11a and lib. In Fig. 11a the target reflectance (mean values^) 
of brittle snow, determined with a reference -pannel, show a good correspondence. 
In contrast to the conformity of the target reflectance is the reflected radiance 
(mean values E) , as shown in Fig. lib with a clear separation between measure- 
ments in sun and in shadow. Similar results were achieved for other snow types. 

The verification tests are currently under way. Since not sufficient data could 
be collected for all main snow types and consequently no definite conclusions 
can be made as yet, the measurements will be continued. 


7.1.2 Mantling of transition zone 

h 

The location of the exact course of the 1 snow line - in., reality a smaller or 
larger transition zone between the totally snowcovered .and snowfrae area - is 
one of the most important problems- in mapping the changes and in measuring the 
areal extent of the snow cover. A detailed study was conducted (STIRREMAHN, 

1977) of a mountainous test site of 2700 km^. (Central Gris op 3 , Fig. 2) on 
this problem, giving special attention to 

* r. :• 

- the effects caused by the relief (snow in sun. and in shadow etc.), 

-/the careful selection ofj the training sampl.es in particular in and around 

the transition zone (Fig. 12), 

- combination with ground data and : aerial photography. 

For an efficient and economic mapping band 5 and 7 were selected and its data 
classified with a PPD-algorithm. 

F j 

The selected categories and. the statistics of the corresponding training samples 
are given in Fig.. 15 and the classification matrix for thes.e training samples 
in Fig. 14. The total area was thgn classified into three m^tin categories; 

- snow oover 

- transition zone. 

- snowfree area. 
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Two types of results were gained: 

- First an exact location of the transition zone (Fig. 15 a and 15 b) in maplike 
form, and easy to compare with an existing topo-map; From this map it is 
possible to extrapolate the exact position of the temporary snow line. 

- Second a calculation of the areal features (Fig. 16): For April 22nd, 1975, 
there were still 65 $ of Central Orisons completely covered with snow and 

9 $ belonged to the transition zone. The accuracy of the classification 
was checked carefully from aerial photography within four sample areas. 

Its findings , su m m a rized in Fig. 17 , conclude that at least 94 ft accuracy 
was achieved. 


7.1.5 System evaluation 

This experiment was carried out by ITTEM (1975) in the US to separate dry and 

n ‘i p 

wet snow, snowcoverei and bare forest with several classifiers. The test area 
was in the Windriver Mtns, Wyoming. The systems used were LARSYS- 3 , STANSORT-II 

* " r 

and IMA.GE-100, The main result out of the comparison is, that human interactions 
throuth direct access to the digital’ operations are an' essential time- and 
money-saving factor. Interactive systems equipped with b&w dr color observation 

'• n i H fi 

screens, to evaluate intermediate results, are to be favourized. It is believed, 

n n 

that a highly interactive, especially designed image processing system (digital 
or hybrid) together with a skilled applications specialist can, for the future, 
bring the best use in satellite snow cover observations. Fig* 18 gives a 
summary of the classification results. 


7,1,4 Outlook and conclusions 


The snow cover project will be developed further in particular in view of 

, _ , . h ^ t A* P O 

change detection to automatically . map the variation of the snow cover as a 
r s*: fc r r 

basic input to runoff models and for water runoff prediction from digital 

:• r. A 

satellite data. This means. additional studies on geometric and radiometric 

, 1 ’5 X 

corrections, on the effects of the fast ohangeing atmosphere, and 'on the cloud 

c .t f J ■ * 

penetration and interpolation problems. It, -must be borne in mind that snow is 

. ' , .1 i j, o *1 

quite a variable feature and subject to strong variations in signature response. 

!< ^ j; ^ t 

As an example the average yalues of all training samples of snow in sun from 
the orbits of Oct, 7th, 1972 and April 22ad, .1975 are compared for the four. 

MSS channels (Figi 19 ), 
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Extension of ground measurements, in particular of the spectral signature pro- 
perties, will lead toward a more accurate classification of the snow cover. 

The addition of the thermal channel on LANLSAT-C could he of great importance 
and interest in this respect. 

Research on the efficiency of the various classification algorithms for 
different purposes will be included, too. Together the analysis will lead to 
the best estimate possible of am. operational monitoring of the areal extent 
of the snow cover. 

A special problem, which occurs in snow mapping stems from saturation of the 
MSS data of fresh snow in sun . Future LANDSAT-HSS sensors should be designed 
to measure snow radiance through the entire range of expected values and not 
be cut off. 

Since clouds will always influence snow mapping in almost all important areas 
of the world a channel in the 1,55 - 1,75 pn band should be considered in 
forthcoming earth resources satellite projects. As learned from ERBP-S-192 
data (MORI, 1976), this channel in combination with one in the green to red 
and one in the near IR region would provide the best -^prerequisites for further 
operational systems (EREP final report, 10. l). 

7.2 Land use mapping for planning purposes 

Various studies. on land use mapping were carried out(and already partly re- 
ported)in different areas and for different purposes. 

7.2.1 Land use and vegetation mapping in the Po Valley 

The results of this investigation (BINZEGGEK, 1975) have been reported in 
NASA Quarterly Report l/2, 1975 (lO.l). 

The main land use. categories considered were: 

- Water 

- Built-up (three. different categories) 

- Agriculture with living vegetation 

- Agriculture with dead vegetation ‘ 

- Forest 
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Por a better classification ratio variables were used in a stepwise, linear 
discriminant analysis* Fig. 20 demonstrates an interesting result of this 
procedure. It is a comparison of two different classifications considering 
only the categories without vegetation (built-up areas, water etc.). The 
ratios of the signals 7/5 are used as a single variable (density slicing), 
where a value smaller than 0,5 denotes vegetation and is presented in blue. 

The classification output for the same categories with discriminant analysis 
is in purple. The combination of the two colors gives the pixels classified 
identically in violett, whilst the blue and purple tones identify the diffe- 
rences. The improvement with the discriminant analysis is obvious: The built- 
up areas and the freeways at the fringe of the City of Milano are delineated 
much better and misclassific&tion in the rural areas are largely eliminated. 

It could be shown that it is possible with LANDSAT-data to either create new 
small scale thematic maps (lilOO’OOO - lsl'OOO'OOO) or to update existing the- 
matic maps. Fig. 21 gives examples of a comparison between the LATOSAT-classi- 
fication and an existing newly editioned topographic map. The detailness on the 
LANDSAT classification is by far superior to the conventional mapping technique. 

7.2.2 Analysis of mult itemu oral data 

To actually map changes in land _use_ & direct comparison of multit amp oral data 
with digital methods is necessary. This asks for additional radiometric and 
geometric corrections. The same principles as discussed in 7.1, will be applied. 

land use studies a test area in' ’the ^Grosses Moos 11 (Cantone of Berne and 
Fribourg, Fig. Z) was selected with -relatively flat topography ; pnd .large 
fields (compared with other parts of Switzerland). 

To superimpose different LAOTJSAT— frames exact reference points have to be 
located. For a precise pixel location mirrors were used to saturate specific 
pixels (4.2). Various circumstances hampered this experiment, in particular 
the problem of receiving the* orbital parameters before the overflight for a 
proper setting of the mirrors (calculation of azimuth and elevation). Thanks 
to the generous support of GSFC we finally succeeded in receiving the data. 

The necessary software was developed to calculate these parameters. The mirrors 

v; 1 1 

could be set up correctly on three cloudfree overflights 'in August and September 
1976, but up to now we have not received the MSS and do fiot know about the 
success or failure of the experiment. 
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Different computer programs are developed to compile the data for different 
overflights by estimating the coordinates of specific pixels on a printer 
output and a corresponding shifting of all pixels (resampling) to achieve 
geometric conformity. Radiometric corrections will be included as discussed 
in 7*1. '4* 

The objectives are to use these corrected multi temporal data 
-to gain additional variables for a supervised classification of the land 
use categories, 

- for change detection in agricultural 4 land use. 


7.2.5 Land use mapping in the Swiss Plateau for planning purposes 

An extensive study is under progress to investigate the possibilities of 

LAUDSAT-data for regional planning purposes. This means not only that the 

results have to meet tough accuracy standards, which have to be achieved 

unter rather difficult conditions (severe relief effects, smallness and heavy 

intermixing of individual features etc.) but also to provide the information 

cheaper than by conventional methods. The prime objective is to delineate as 

exactly as possible the extent of the settlements and its separation from 

the surrounding open land. The test area includes the total Swiss Plateau, 

e.g. the heavily populated area of Switzerland and covers approx. 10*000 km (Fig. 2). 

Processing is done with the IBIS system (Fig. 6) using a PPD classifier. 

To achieve this goal numerous tests have Been executed lately, partly by 

t 

< students under the leadership of Dr.K. Itteh. This includes; 

- Land use studies around the Lake of Pfaffikon (Kt. Zurich). Here the problems 

, 1 ' fi n 

are to achieve a good classification of smaller settlements and the various 
rural land use categories. The separation of open water surfaces, wetlands 
and meadow/paatures were given specific attention. As it turned out it was 

not possible on a fall image to separate wetlands from settlements. The 

" ^ ,, , 

"city" in the middle of the wetlands at the Southern shore of Lake Pfaffikon 
is a miscla3sification, which could not be eliminated (Fig. 22). 

- Analysis of different urban areas of various size by 

• 5 » v 

= a careful study of the fringe zone 

- differentiations within the settlement into high density built-up areas 
and low density built-Up areas. • 

These studies were carried out for the Cities of Schaffhausen, Winterthur, 
Frauenfeld, Rapperswii, Zug, .Baar, Cham,- Luzern (Pig. 2)y with a population 



- 21 - 


between 110*000 and 8 * 000 . In Fig. 23 the three categories "water, 
high density built-up and low density built-up" are presented for the 
City of Schaffhausen. 

- Land use of the agglomeration of Luzern (Fig. 24) and Zurich (Fig. 25 ). A 
general classification was undertaken with special reference to the above- 
mentioned problems again around cities of various size.’ Eight different 
categories were distinguished. 

Sofar it may be concluded that a fall imagery (October 7th) is not suitable to 
provide a clear separation of built-up areas and wetlands, • since the spectral 
range of these two categories is almost identical. (Therefore a spring picture 
has to be used in addition (test will be undertaken on a frame of April 22 nd, 
1975) and combined with the fall image. 

After the completion of these tests and a careful evaluation of the training 
samples the total study area will be classified and the results presented with 
the photomation system. 

7.2.4 Correlation between size of settlement and population 

An interesting study on this topic was undertaken by a group of ' geography 
students of the University of Zurich (GALLUSER et al, 1976) for the Swiss 
Plateau (e.g. the heavy populated part of Switzerland) using two different 
methods. * 

The first method, consists of delineating and measuring the areal extent of all 
identifiable settlements. For the Swiss Plateau official planning organizations 
have evaluated that the average value for closed settlements is 100 persons per 
ha. Since streets, squares and other parts without buildings within a city 
cannot be eliminated, the average measurement of the areal extent is about 30 $> 
too high. But this value compensates all small and individual settlements, 
which cannot be identified in a satellite image. 

• i * 

Tests were undertaken for the agglomeration of 

- Basel (Fig. 26): The estimation was 4,2 $ too high. 

2 , 

- ¥interthur: The estimation was 5 % too small. 

- Four small rural settlements: The estimation was 4,68 $ too high. 
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Afterwards the test areas, Kt. Zurich and Kt. Schaffhausen (Fig. 2) , were 
interpreted with special reference to the three big agglomerations. 

The second method is based on a delineation of all forests and water bodies. 

The remai nin g area (BRUTTCFLAECHS - BF) can be measured and correlated with 
existing population values from planning organizations and statistical year- 
books. The corresponding values for the test area are: 

- Kt. Zurich: 10 Persons / ha BF 

- Kt. Schaffhausen: 4,5 Persons / ha BF 

With this method the settlement itself has not to be recognized and delineated, 
•the only problem is to not mix it with forests. The method consequently is 
very suitable for areas with small and/or scattered settlements. 

The following results were' achieved (average of results from three students) 
from satellite data of August 9th, 1975; the population statistic dates from 
1972. 


/ r» 


First method 
Second method 

Population estimation of Kt .Zurich and Schaffhausen 

.+ 1,6 

• . - 12 % 



' First method 

Agglomerations (statistical data 1974-)' 

Zurich 
- 7,3 i> 

Winterthur 
' + 15,5 % 

A y 

Schaffhausen 
+ 7,2 £ 


The results show a surprising good estimation for the first method, whilst 
for the second one an unexpectedly big error occurs. Of the Jhree agglome- 
rations Winterthur brought the biggest divergence, which is.. easy to explain 
by the rugged topography with large forests in the immediate’ surroundings 
of the city. This handicaps an exact delineation of the built-up areas. 
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X».2,5 Outlook and conclusions 

Sofar the possibilities of classifying land use features and the creation of 
SI&&2.X scale h em c rnaips could, Iso demons i; rated successfully* 

Similar problems as discussed for snow mapping occur when making the systems 
oruely operational for change detection. It can be achieved by extending 
additional geometric and radiometric corrections'.', 

The exact delineation of the fringe zone of settlements of different sizes 
under European condition is still subject of further investigations. 

Correlations of areal extent of settlements and standard population indices 
from planning agencies provide good results. 

Satellite data will be of great interest to. regional planning purposes , even 
in a small and heterogeneous country like Switzerland. Research work there- 
fore is directed toward these applications to provide planning agencies 
with the latest situation on the ground,. in particular of recent changes 
m the built-up areas, and in correspondence with movements, of the popu- 
lation. 



- 24 - 


8. .FUTURE ACTIVITIES 

The investigation will he continued in the following directions, combining 

the three aspects 

- application of satellite data for hydrological (snow mapping) and planning 
purposes (land use studies) , 

- data processing, 

- interactive data analysis. 

a) Snow mapping of large areas in the Alps considering extensive radiometric 
and geometric corrections. 

b) Snow mapping in well defined watersheds. The areal measurements will be 
used as input for an existing runoff model. 

c) Modification of selection and evaluation of training samples using an ' 
interactive system. 

d) Increased correlation of data from ground' (EXOTECH measurements) and air- 
borne scanners (l^S) with satellite data (atmospheric measurements, corre- 
lation with models of standard atmospheres etc.). 

i 

e) Different at ion. of the snow surface (snow type, surface conidi tions, 

0° C temperature line dtc . ) . 

f) Change detection in rural areas. Digital processing of multitemporal MSS- 
data for land use studies,' considering geometric corrections and resamp- 
ling techniques. 

g) Improvement of' soft ware 0 systems' 1 regarding 

- reformating 'of data (transformation of rectangular pixel format into 
a square form), 

= correlation of multrcemporai- data, 

= change detection. 

h) Improvement of control programs for photomation system (transfer functions 
etc. ) . 

i) Combination of] aatellitp data with existing data from other sources (geo- 
graphical information systems)..* 

k) Evaluation of ’various algorithms for specific topics (cost/benefit studies). 



l) Application of satellite data for regional planning, e.g. changes in 
rural and urban land use and connected migration of population. 

m) " Analysis of settlements of various sizes and functions especially regar- 

ding a detailed delineation of its fringe zone. 

For various of the above mentioned problems, especially on snow studies, 
the availability of the LA.NDSAT-C thermal scanner data will be of parti- 
cular interest. 



9* SIGNIFICANT EESULTS 


A system was developed for operational snow and land use mapping, based on 
a supervised classification method using various classification algorithms 
and representation of the results in maplike form on color film with a 
photomation system. 

In particul it could be shown that snow mapping in high mountain terrain 
can be achieved independently of relief and shadow effects, sun angle, snow 
type, surface conditions etc. The transition zone between the totally snow- 
covered and the totally snowfree areas, which is of specific interest in the 
melting process and for runoff predictions, was analysed in detail regarding 
location, extent and vertical dimensions, with an accuracy of approx. 94 ft. 
To separate snow from clouds and .the snowfree background the four IANB3AT- 
bands are insufficient and should l be supplemented in future by a band in , 
the 1,55 - 1,75 pm regions(as learned from ERSP-S-192 data}* A separation 
of various snow types, especially of fresh and old snow', will be possible as 
learned from BXOTECH measurements on the ground* The ratios 4/7 end also 
5/7 are good indicators for the snow moisture conditions. 

Land use mapping under European, conditions (relief effects, small size of 
individual, fields^ phenologieal differentiation, cultivation 1 " practices, 
heavy intermixing ' of many drop types etc.) could he achieved with a stepwise 
linear discriminant analysis by using additional ratio >var±ables. For regio- 
nal planning purposes special - attention was given to a 'detailed delineation 
> ^ 

of the built-up areas. On 'fall images these signatures i are often hot sepa- 
rable from wetlands. Two different methods were tested i.to correlate the 

* 

size of settlements and the population with ah accuracy) for tie densly popu- 
lated Swiss Plateau betweeh + 2 and - 12 ■*> 
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LAl'TDSAT coverage of Switzerland 1975 - 1976 as received from 
GSFC and Telespazio, Rome 
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Fig. 1; (continued 1 ) 
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Fig. 3: Example of high altitude aerial underflight; Davos and vicinity, 

November 24th, 1975 

(photo: Swiss Air Reconnaissance, DUbendorf) 
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Fig. 4: Decision making techniques in digital processing of multivariable 

data (after STSIUSE, 1972) 
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IBIS interactive interpretation system for IBM 370/155 (after FASLER, 1976) 
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Pig. 7 : Reflectance g {%) of various snow types in the 4 LAND3AT MSS 

channels measured with EX0TECH-100 (-after STAENZ, 1976) 



CHANNEL 


. Legend PA = old powdery snow 
P = sticky snow 
B = slush snow 


BR = brittle snow 
S = spring snow 
TK = solid crust 
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‘ig. 8 + 9: Target reflectance ($>) as function of snow depth 3H (cm) measured 

with EXOTECH- 100 (after STAENZ, 1976) 

Pi g. 8: channel 4 
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Fig. 10: Ratio values (RW) of target reflectance of EXOTECH-measurements . 

in function of time during melting of snow; near Davos, Grisons, 
May 20th, 1975 (after STAENZ, 1976) 
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Fig. 11: Comparison of reflectance measurements in sun and shadow 

(after STAENZ, 1976) 



left: target reflectance ( ^ ) of brittle snow determined with reference 
panne 1 

right: Radiance (I) 

A = in sun 
B =s in shadow 



Fig. 12: 
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Selection of training sample for snow near transition zone 
(after STIRNEMANN, 1977) 
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Statistics of training groups {mean values and standard deviation) 
(after STIEHBMAHK, 1977) 
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Classification of -training groups, snow mapping. Central Grisons, April 22nd,- 1975 
(after STIRNBI-IANN, 1977) 
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Fig. 15 a+b : Digital snow mapping of Central Grisons for April 22nd, 1975 

(after STIRH3&1ANN, 1977) 



a) Map of total test area 
white = snow cover 
gray = transition zone 
black = snowfree area 


beproducbtlek of 

ORIGIN All PAGE IS POOB. 



b) Enlarged section showing 
location and extent of 
transition zone • 
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Fig. 16: Classification and areal measurement of total test site 

(after STIRWEMM, 1977) 
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Fig* 17 1 Accuracy of LANDSAT-classification compared with air photo 
measurement for 4 test 3ites (after STIRNEMANN, 1977) 
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Fig. 18: Comparison of results received from 3 different image interpre- 

tation systems (after ITTEN, 1975) 
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Fig. 19 : Comparison of video value of snow in sun from different LA.NDSAT- 

orbits (mean values of training samples in the 4 MSS channels) 
(after STIRNEMAM, 1977) 
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Fig. 20: Comparison of land use classification with density slicing and 

linear discriminant analysis using ratio variables 
(after BINZiSGGER, 1975) 



purple: built-up areas classified with discriminant analysis 

blue: with density slicing 

violett: identical classification 


asssffssws- 
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fig, 21; Comparison of LAHDSAT-classification with new topographic maps of 
similar scale (after BINZEGGER, 1975) 
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Fig. 22: Land use mapping around Lake of Ffaffikon, October 7th, 1972 

(after ITTEN, 1976) 
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Delineation of high-density and low-density built-up areas 
versus open land. Schaffhausen and vicinity, October 7th, 
1972 (after ITTEN, 1976) 
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Land use mapping of the agglomeration of Lucerne, October 7th, 
1972 (after ITTEN, 1976) 
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REPRODUCIBILITY OP THE 
ORIGINAL PAGE IS POOR 


Fig, 25: Land use mapping of the agglomeration of Zurich, October 7th, 

1972 (after* ITTEK, 1976) 
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Pig. 26 a+b: Correlation between size of settlement and population: 

Agglomeration of Basel, August 9th., 1975 (after GALhlJSSER 
et al, 1976) 


a) 


Measurement of the residential area on LARDSAT MSS 7 
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b) Results from areal measurements from LARDSAT and Topographic Map l:25 f 000 
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Diff. 
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Topo , Map 
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Diff. 
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area 

3' 975 ha 


3 '710 ha 


population 

397' 500 

+ 4>2 % 

371' 000 

- 2,7 

Statistic 

1970 

-381 '453 





